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ABSTRACT 

The weak, R-type ionization fronts (I-fronts) which swept across the intergalactic 
medium (IGM) during the reionization of the universe often found their paths blocked 
by cosmological minihaloes (haloes with virial temperatures Tvir ^ 10^ K). When this 
happened, the neutral gas which filled each minihalo was photoevaporated; as the 
I-front burned its way through the halo, decelerating from R-type to D-type, all the 
halo gas was eventually blown back into the IGM as an ionized, supersonic wind. In 
a previous paper (Shapiro, Ihev & Raga 2004, hereafter Paper I), we described this 
process and presented our results of the first simulations of it by numerical gas dy- 
namics with radiation transport in detail. For illustration we focused on the particular 
case of a IO^Mq minihalo which is overrun at z = 9 by an intergalactic I-front caused 
by a distant source of ionizing radiation, for different types of source spectra (either 
stellar from massive Pop. II or III stars, or QSO-like) and a flux level typical of that ex- 
pected during reionization. In a Cold Dark Matter (CDM) universe, minihaloes formed 
in abundance before and during reionization and, thus, their photoevaporation is an 
important, possibly dominant, feature of reionization, which slowed it down and cost it 
many ionizing photons. In view of the importance of minihalo photoevaporation, both 
as a feedback mechanism on the minihaloes and as an effect on cosmic reionization, 
we have now performed a larger set of high-resolution simulations to determine and 
quantify the dependence of minihalo photoevaporation times and photon consumption 
rates on halo mass, redshift, ionizing flux level and spectrum. We use these results to 
derive simple expressions for the dependence of the evaporation time and photon con- 
sumption rate on these halo and external flux parameters which can be conveniently 
applied to estimate the effects of minihaloes on the global reionization process in both 
semi-analytical calculations and larger-scale, lower-resolution numerical simulations 
which cannot adequately resolve the minihaloes and their photoevaporation. We find 
that the average number of ionizing photons each minihalo atom absorbs during its 
photoevaporation is typically in the range 2-10. For the collapsed fraction in mini- 
haloes expected during reionization, this can add « 1 photon per total atom to the 
requirements for completing reionization, potentially doubling the minimum number 
of photons required to reionize the universe. 

Key words: hydrodynamics — radiative transfer — galaxies: halos — galaxies: high- 
redshift — intergalactic medium — cosmology: theory 



1 INTRODUCTION 

The Cold Dark Matter model, due to its hierarchical na- 
ture, predicts that a significant fraction of the matter in the 
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universe before and during reionization (up to 30 — 40% 
at redshift z — 6 for current ACDM) resided in collapsed 
and virialized minihaloes - small haloes with virial tempera- 
tures Tvir below 10* K. These minihaloes, in the mass range 
between the Jeans mass (below which gas pressure forces 
prevented baryonic matter in the IGM from collapsing into 
dark matter haloes as they formed; i.e. Mj < W^Mq) and 
the mass for which Tvir = 10* K [M(10* K) < 10** Mq] were 
filled with neutral gas, since their temperatures were too low 
for collisional ionization to be effective. Cosmic reionization 
had a profound effect on the gas content of these minihaloes 
and they, in turn, interfered with cosmic reionization. 

This reionization began when the first sources of ioniz- 
ing radiation turned on, causing global I-fronts to sweep out- 
ward through the surrounding IGM, creating intergalactic 
H II regions, and ended when these H II regions became large 
enou gh and numerous enough to overlap iShapiro fc Giroua 
Il987t) . These global I-fronts were generally weak, R-type 
fronts, which moved supersonically, not only with respect to 
the cold, neutral gas ahead of them but also with respect to 
the gas photoheated to T>10*K behind them. As such, they 
raced ahead of the hydrodynamical response of the ionized 
gas, unaffected by the disturbance they created in the gas be- 
hind them. That situation changed, however, when the path 
of one of these global I-fronts was blocked by a minihalo. 
Since the speed of an I-front varies inversely with the density 
of the neutral gas ahead of it, the I-front must have slowed 
down when it encountered the dense, centrally-concentrated 
gas inside a minihalo. The latter is, on average, more than a 
hundred times denser than the mean IGM, while at its cen- 
ter it is yet another hundred times denser even than that. 
This denser gas slowed the I-front down enough to "trap" 
it, long enough for the gas-dynamical back-reaction of the 
disturbed minihalo gas to catch-up to the I-front and affect 
its progress. Thereafter, the fate of this I-front and that of 
the minihalo gas were inextricably coupled, as ionized mini- 
halo gas, heated to T ^ lO^K, above the binding energy of 
minihalo gravity, blew back toward the ionizing source, into 
the IGM as a supersonic, evaporative wind. 

The first realistic discussion of the photoevaporation 
of cosmological minihaloes by global I-fronts during reion- 
ization, including the first radi ation-hydrodynamical sim- 
ulations of this process, was by Sha^irOj^^a ga^ Mellemal 
Jl997lll998ri. de v eloped further in lShaoiro fc Raeal ll2000allbl 
I2OOII) : IShapird J200ll) and Paper I (see Paper I for a 
more complete description and additional references.) As 
described in Paper I, the dominant sources of reionizing 
radiation are likely to have been haloes of virial temper- 
ature Tvir ^ 10*K. Star formation in minihaloes, instead, 
is thought to have been inhibited early on by the pho- 
todissociating effect on minihalo H2 molecules of the ris- 
ing UV background just below the H Lyman continuum 
edge, created by a small fraction of the sources ultimately 
required to complete reionization. If so, then each source 
of reionization would typically have found its sky cov- 
ered by minihaloes between it and the nearest neighbor- 
ing source haloes, and the expansion of I-fronts in the 



IGM around each source to the point of overlap would 
have depen ded upon the photoevaporation of intervening 
minihaloes (Haiman, Abel & Madau "2001"; 'Shapiro' '20011; 
[Barkana fc Locb 2002; Shapiro. Iliev. Raga & Martel 2003; 
IPapeiMl Tsince the gas inside collapsed minihaloes was much 
denser than average, the rate of recombination there was 
much higher than average, as well. Hence, minihalo pho- 
toevaporation was likely to have consumed multiple ioniz- 
ing photons per atom, thereby increasing, and potentially 
domina ting, the global photon consumption during reion- 
ization jHaiman. Abel fc Madaull200ll : IShapirolEoOll : IPaperl 
Q). Therefore, quantifying this process is crucial for deter- 
mining the onset, progress and duration of the reionization 
of the universe. 

An early start (2:>15) and late finish (2 ~ 6) for cosmic 
reionization are suggested by recent observations of the fluc- 
tuating polarization of the CMB ( Kogut et al. 2003) and 
the Gunn-Peterson effect in the spectra of quasars at z ~ 6 
(Becker et al. 2001, Fan et al. 2003), respectively. This is 
consistent with the view that the volume-filling factor of ion- 
ized regions of the universe grew over time from 15 until 
they overlapped finally at z « 6. The presence of minihaloes 
filled with neutral gas was affected by this reionization, not 
only because photoevaporation stripped pre-existing mini- 
haloes of their baryons, but also because the subsequent 
collapse of baryons into newly forming minihaloes was sup- 
pressed by "Jeans-mass filtering" of the linear density fluctu- 
ations in the baryonic component after the IGM was heated 
to T > W^K when it was reionized (Shapiro, Giroux, and 
Babul 1994; Shapiro 1995; Gnedin 2000b). These two pro- 
cesses acted to reduce the presence of neutral-gas-filled mini- 
haloes only in the ionized regions, however. The neutral vol- 
ume outside these regions continued to form new minihaloes 
without interruption, at the unflltered rate of the universe 
without reionization. As a result, the I-fronts which led the 
expansion of the ionized regions would have advanced into 
neutral gas in which minihaloes formed with ever- increasing 
abundance over time. The unflltered collapsed fraction in 
minihaloes for ACDM ai z — 20, 15, 9, and 6 was 3, 10, 
28, and 36 percent, respectively (Paper I). The question of 
whether minihalo photoevaporation was a dominant feature 
of reionization will ultimately be answered only when we 
know the full story of cosmic reionization, since there also 
exist (somewhat more speculative) scenarios in which the 
formation of some minihaloes might have been suppressed 
by, for example, a signiflcant X-ray background at high-z 
(e.g. Machacek, Bryan & Abel 2003, Glover & Brand 2003, 
Mad au et al. 2004, Ricotti & Os triker 2004), double r eioniza- 
tion (ICenl i2003). entropy floor JOh fc Haimanll2003^ . or re- 
duced small-scale po wer in the power spectrum of primor dial 
density fluctuations llSomerville. Bullock, fc Livioll20o3l . 

In Paper I, we discussed the process of photoevapora- 
tion of cosmological minihaloes in detail. We demonstrated 
the phenomenon of I-front trapping, whereby the highly su- 
personic, weak, R-type front propagating in the low-density 
IGM slows down upon entering the minihalo, where the gas 
density is much higher than that of the IGM, and converts to 



© 0000 RAS, MNRAS 000, 000-000 



Photoevaporation of cosmological minihaloes II 3 



D-critical and then subcritical D-type, preceded by a shock 
which compresses the gas ahead of the I-front. The I- front 
slowly propagates subsonically through the halo, stripping it 
of gas, layer by layer, as this gas, once ionized, expands into 
the IGM as a supersonic wind towards the ionizing source. 
This process eventually results in a dark matter halo com- 
pletely devoid of gas. We described there in detail the struc- 
ture and evolution of the photoevaporative flows and their 
dependence on the radiation spectrum of the external ioniz- 
ing source, considering both massive stars (Pop II and Pop 
III) and QSO-like spectra. We derived the evolution of the 
position and speed of the I-front, the evolution of the halo's 
remaining neutral mass fraction and effective opaque geo- 
metric cross-section for the absorption of ionizing photons, 
as well as the global parameters of photoevaporation like 
the photoevaporation time, t^v, and the number of ionizing 
photons absorbed per atom in the course of that evapora- 
tion, 1^, and its evolution. Finally, we also discussed some 
observational diagnostics of the photoevaporation process. 

In this paper, we extend the results of Paper I by in- 
vestigating the dependence of the photoevaporation process 
on the mass of the halo, on the level and spectrum of the 
external photoionizing flux responsible for the I-front which 
overtakes the halo, and on the redshift at which this I-front 
encounters the minihalo. As discussed in Paper I, we make 
the assumption that the H2 molecule formation and cooling 
which might have led to star formation inside minihaloes 
were inefficient in the bulk of the minihalo population, so 
the minihaloes were largely "sterile" reservoirs of neutral 
atomic gas. We concentrate particularly on the dependence 
of the two properties of minihalo photoevaporation which 
are most important for quantifying their global effect on 
cosmic reionization, namely, the evaporation time tcv and 
the net ionizing photon consumption rate per minihalo (i.e. 
the number of ionizing photons absorbed per minihalo atom 
during the evaporation time), ^. In §|21 we shall define the 
quantities tcv and ^ and how we use our simulations to eval- 
uate them more precisely, and recall some approximations 
which have been used previously to estimate these quan- 
tities, for comparison. Our numerical simulations are de- 
scribed briefiy in §|3 In S 14. II and S 14.21 we summarize our 
numerical results for tcv and ^, respectively, for haloes of 
different masses encountered at different redshifts by the I- 
fronts driven by external ionizing sources of different flux 
levels and spectra, and compare these with the previous es- 
timates described in § |5| In § 14.31 we show that the time- 
dependence of the evolving neutral mass fraction inside a 
minihalo during its evaporation has a universal shape, in- 
dependent of halo parameters and external flux level. The 
implications of these results for the theory of cosmic reion- 
ization are briefly discussed in §|S] Throughout this paper we 
use the current concordance model for the background uni- 
verse and the dark matter - flat, COBE-normalized ACDM 
with Slo = 1 - Ao = 0.3, h = 0.7 and fib^^ = 0.02, with 
primordial density fluctuation power spectrum P(k) oc fc"'', 
with Up = 1, and Ugi^-i — 0.87 (e.g. Spergel et al. 2003). 



2 HALO EVAPORATION TIMES AND 
IONIZING PHOTON CONSUMPTION 
RATES 

We define the evaporation time tov to be the time at which 
only 0.1 per cent of the gas which was initially inside the 
minihalo remains neutral^. As in Paper I, we then obtain 
the number of ionizing photons absorbed per minihalo atom 
during this evaporation time by directly counting the num- 
ber of recombinations experienced by each of those initial 
halo atoms and its first ionization: 

m = ^ + 7^ / '^^("h^'^hii + a'ilnnou)n,, (1) 

where Ue is the number density of electrons, nnii and riHeii 
are the number densities of H II and He II, respectively, while 

and a^^' are the Case B recombination coefficients for 
H II and He II, respectively, A^a is the total number of atoms 
initially inside the minihalo, and Mon is the number of these 
atoms initially inside the minihalo which are ionized by the 
evaporation process. At time t = tev, practically all atoms 
originally in the minihalo are ionized, thus A'^ion = Na and 
^ = C(iov) is the total number of ionized photons consumed 
per minihalo atom. We have neglected the recombinations 
of He HI to He II because these generally contribute diffuse 
flux whi ch is absorbed on the spot by H and He. 

Hai man. Abel fc Ma dau (2001) estimated the evapora- 
tion times and photon consumption rates for evaporating 
minihaloes analytically as follows. The evaporation time tcv 
was approximated by the sound-crossing time for the char- 
acteristic size of the minihalo at the sound speed of the ion- 
ized gas, tsc ~ 2rt/cs{W'^K), roughly consistent with our 
earlier results for th e simulation of the photoevapora t ion of 
a lO'^Mrr, minihalo ishapiro. Raga fc Mellernal Il99l Il998l: 
IShapiro fc Rag3l2000allbll200ll : IShapirdl200 J) . Haloes were 
assumed to be dark-matter dominated, nonsingular isother- 
mal spheres a ccording to the T runca ted Isothermal Sphere 
(TIS) model JShapiro. Ihev fc RagalHoQi llliev fc Shapird 
I2OOII) . For a minihalo of mass M at high redshift z, the 
TIS halo solution yields 

t.. = 98Myr(M.)^/^ (tT^) + (2) 

using the adiabatic sound speed Cs{T) — 11. 7 (T4, / fi)^^^ — 
15.2Tl^\ms-^, where fj, = 0.59 is the mean molecular 
weight for fully ionized gas of H and He if the abundance 
of He is j4(He) = 0.08 by number relative to H. Based on 
this evaporation time estimate, a naive estimate of ^ is then 
made possible by further assuming that the minihalo has 
zero optical depth to ionizing photons and is fully ionized 
instantaneously but remains static at its initial density for a 
time tsc (the optically-thin, static approximation, hereafter 



^ As discussed in Paper I, the evaporation time is not sensitive 
to the precise value adopted for this remaining neutral fraction 
in the definition, as long as the latter is much less than unity. 
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referred to as "OTS"). Ignoring the contribution of He to 
recombinations, equation Q then yields 



3 THE CALCULATION 



. 1 , , Cint{nH)a^H\ 

^OTS = 1 + / ——z 

i + Otis 



(3) 



jHaiman. Abel fc Madal 1200 ill (where we have added the 
first term on the r.h.s. above to their equation to account 
explicitly for the fact that each atom must first be ionized 
once before it can recombine and be ionized again), where 
{tih) is the mean H atom number density inside a halo, 
Cint = (nil) / {nn}'^ = 444^ is the effective clumping factor 
for the TIS, and 1 + iStis = 130.6 is the average overdensity 
of a TIS halo with respect to the cosmic mean background 
density. The quantity / is an "efficiency factor" introduced 
to take account of details neglected by the OTS assumptions. 
Together, equations ||5J and @ yield 



Cots = 1 + 206 f T-^^* M^^"" 



0.15 



-1/3 



(4) 



Here, we have defined il/7 = Af/lO^A/o and (1 + 2)10 = 
(1 + z)/10. ^aiman. Abel & Madau (20Q1,) cahbrated the 
factor / against an optically-thin simulation of the expan- 
sion of a uniformly photoheated halo (i.e. no radiative trans- 
fer), obtaining / ~ 1. According to this OTS approximation, 
Cots S> 1 for minihaloes during cosmic reionization, thus 
the presence of minihaloes during cosmic reionization would 
have dramatically increased the nu mber of ionizing photons 
neede d to complete reionization iHaiman. Abel fc Madaul 
120011) . According to equation that is, each atom which 
had collapsed into a minihalo before reionization ended 
would have consumed many ionizing photons during its pho- 
toevaporation, a disproportionately large share of the ion- 
izing photon background compared to that which would be 
available to ionize the uncoUapsed atoms in the diffuse IGM. 

As shown by the detailed simulations in Paper I, for 
haloes of mass lO^M© exposed at z; = 9 to the typical level 
of ionizing radiation expected during reionization, while the 
OTS approximation yields a reasonable, rough estimate of 
tev, it grossly overestimates ^, by factors of 30-50, depend- 
ing on the ionizing source spectrum. This more accurate 
determination of C in Paper I by numerical hydrodynami- 
cal simulations with radiative transfer nevertheless showed 
that minihalo photoevaporation could still have consumed 
enough photons to increase substantially the total number 
required per baryon to complete reionization. By contrast, 
earlier estimates of the requirements for reionization of the 
IGM had concluded that only about one photon per atom 
was required, but these estimates were based on calculations 
without sufficient resolution to account for the m inihaloes 
JCnedinlbOOOatlMiralda-Escude. Haehnelt fc Reej EoOO). It 
is therefore crucial to extend the accurate determination of 
C in Paper I to the full range of minihalo masses, ionizing 
flux levels and redshifts expected during reionization, in or- 
der to evaluate the global effect of minihaloes on the process 
of reionization. 



3.1 Numerical Method and Initial Conditions 

We have performed a large number of simulations of the pho- 
toevaporation of individual minihaloes which occurs when 
a weak, R-type intergalactic I-front encounters a minihalo 
during cosmic reionization. Our simulation method and ini- 
tial conditions were described in some detail in Paper I, so 
we only briefly summarize these here. Each halo is modelled 
as a nonsingular TIS (Iliev fc Shaoiro 2001), surrounded by 
self-similar c osmological infall according to the solution of 
iBertschingeJ lIlQSgf) . We generalized this self-similar infall 
solution to apply also to a low-density ACDM b ackground 
universe at high redshift jlliev fc Shapiro! I200I . At high- 
z, the halo radius for the current background cosmology 
adopted here is rt — 0.754 M^''"^ (I -I- Zcoii)]~o^ kpc, the to- 
tal atomic number density at the center is no = 3.2(1 -|- 
-ZcoiOio cm""^, and the halo virial temperature and velocity 
dispersion are Tvir ~ 4000M7^^(1 -I- Zcoii)io K and av ~ 
5.2A'/7''^(1 + 2coii)io^ kms~^, respectively. For haloes with 
Tvir ^ lO'^K (minihaloes), the gas atoms are neutral. We 
consider haloes in the mass range between the Jeans mass 
of the uncoUapsed IGM prior to reionization, Mj = 5.7 x 
10' (f7o/iV0.15)"'^' X (n6/iV0.02)"'''' [(1-f z)^o]'''' AfQ, 
and the mass for which Tvir ~ IC^K according to the TIS 
model, A'/(10''K) = 3.95x 10^(J7ohV0.15)~'/^[(l-hz)io]-'/^ 

Our simulations use the two-dimensional, axisymmet- 
ric Eulerian gas-dynamics code CORAL, with radiative- 
transfer. Adaptive Mesh Refinement (AMR), and the van 
Leer Flux- Vector Splitting Algorithm, which resolves shocks 
well and properly tracks I-fronts ranging from fast, super- 
sonic R-type to sub-critical, subsonic D-type, as described 
in detail in Paper I (and references therein). All our simu- 
lations were set up as we described in detail in Paper I. We 
fixed the size of the computational box in all simulations 
presented here to be the same in units of the halo radius 
(with the long side ~ 10 times larger than the halo radius), 
Xbox oc Tvir cx: M,|{j^(l -|- Zcoii)"^, SO that the fully-refined 
grid resolution is equivalent in all cases. As in Paper I, we 
consider three different cases for the spectrum of the exter- 
nal source of ionizing radiation: (1) starlight with a 50,000K 
black-body spectrum (hereafter, referred to as "BB5e4"), 
representative of massive Population II stars; (2) starlight 
with a 100,000K black-body spectrum (hereafter, "BBle5"), 
as expected for massive Population HI stars; and (3) QSO- 
like, with a power-law spectrum with spectral index —1.8 
(hereafter, "QSO"). As in Paper I, we express the unatten- 
uated external flux of ionizing photons as a dimensionless 
quantity, Fq, the flux in units of that from a source emit- 
ting TVph = lO^^s"^ ionizing photons per second at a proper 
distance d of 1 Mpc, or equivalently, emitting lO^'^ s~^ at a 
distance of 10 Kpc, so _Fo 

^Mpc- 



F/{10''^s-V[47r(lMpc)2]} 
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Table 1 . Numerical convergence test for the evaporation time tev and photon consumption rate per atom § as functions of grid resolution 
for minihalo of initial mass 10^ M©, exposed to initial flux of _Fo = Ij starting at z = 9. 



spectrum 


N 11 


(Ax) ■ foci 


' t / V /mm 


t„„ fMvrl 


f 
s 


BB5e4 


128 X 256 


27.0 


23.4 


465 


6.41 


BB5e4 


256 X 512 


13.5 


56.7 


300 


5.73 


BB5e4 


512 X 1024 


6.7 


113 


195 


5.39 


BB5e4 


1024 X 2048 


3.4 


227 


150 


5.14 


QSO 


128 X 256 


27.0 


23.4 


112 


6.20 


QSO 


256 X 512 


13.5 


56.7 


103 


5.45 


QSO 


512 X 1024 


6.7 


113 


100 


5.20 


QSO 


1024 X 2048 


3.4 


227 


97 


4.95 


BBle5 


128 X 256 


27.0 


23.4 


215 


5.16 


BBle5 


256 X 512 


13.5 


56.7 


178 


4.03 


BBle5 


512 X 1024 


6.7 


113 


145 


3.39 


BBle5 


1024 X 2048 


3.4 


227 


128 


3.33 



3.2 Numerical Resolution 

To begin, we have studied the numerical convergence of our 
simulations with increasing spatial grid resolution, for the 
quantities tev and ^, in order to establish the limits of ro- 
bustness and reliability of our results. We have performed a 
series of four simulations, with fixed halo parameters and 
external flux level, with our fully-refined grid resolution 
NccUs = Nr X Nx, in 2D, axisymmetric (r, a;)-coordinates, 
increasing by factors of 2 from 128 x 256 up to 1024 x 2048, 
for each of the three incident spectra. The results are shown 
in Table0 We note that with increasing resolution both the 
evaporation time t^v and the number ^ of ionizing photons 
consumed per atom decrease, leading ultimately to conver- 
gence in both quantities at our highest resolution. The con- 
vergence is more readily achieved for the cases with harder 
spectra (QSO and BBle5) than for the softer spectrum case 
BB5e4. In all cases ^ converges faster than tev 

3.3 Parameter space 

In order to study how the photoevaporation of minihaloes 
by global I-fronts during reionization depends upon the 
halo mass, the external ionizing fiux level and spectrum, 
and the redshift at which the I-front encounters the mini- 
halo, we have performed the following large set of simu- 
lations. We considered a range of halo masses from a low 
of 1O*M0, which is close to the Jeans mass at the corre- 
sponding epoch, to 4 X 1O^M0, close to M(10^ K), the mass 
at which the halo virial temperature is T^ir ~ 10'' K (e.g. 
M(10*K) = 4 X lO^Mo at z=9). The intergalactic I-front 
which encountered the minihalo in each case was assumed 
to reach the boundary of our simulation volume at a start- 
ing redshift Zi (and we assumed Zcoii ~ Zi in defining the 
initial parameters of the TIS model for the halo). We took 
Zi to range from 1 + Zi = 20, consistent with the early start 
for reionization implied by the recent WMAP results, to 
1 -|- «i = 7, at which epoch the data from QSO absorption 
line spectra suggest that reionization was just ending. The 
flux levels assumed ranged from Fq — 0.01 to Fq — 10"^. 



We argued in Paper I that the interval Fq = 0.01 — 100 
corresponds approximately to the typical range of fluxes ex- 
pected when a global I-front encountered minihaloes dur- 
ing the reionization epoch. Here we extend this range up to 
Fo = 10^, to accommodate cases of very high flux due either 
to a rare, very strong source or to the minihalo's proximity 
to the source. In all our simulations the ratio of box size 
to halo radius was kept fixed at I/box ~ 10, which proved 
sufficient to follow the evaporation fiow until full evaporation 
was achieved. Such a box size guaranteed that our transmis- 
sive boundary conditions were self-consistent at all times, 
since the simulation volume (fixed over time in proper co- 
ordinates) was always large enough to contain the sphere 
of radius equal to the turn-around radius in the cosmolog- 
ical infall profile centered on the minihalo as that radius 
increased over time. We show the complete set of simulation 
parameters and results for tcv and ^ in Figure In order to 
make such a large set of simulations feasible in terms of both 
computing time and storage requirements, we limited the de- 
gree of refinement of the mesh so that the fully-refined mesh 
had 512 x 1024 cells in all cases. Our numerical convergence 
tests above showed that, at this resolution, the simulations 
are already largely converged. To correct for any remaining 
small differences from the fully-resolved limit, we used the 
highest-resolution simulations in Table0to re-normalize our 
plotted results for t^v^ . Furthermore, in order to verify this 
approach, we performed a limited set of additional higher- 
resolution (1024x2048 fully refined) simulations for different 
halo masses and external fluxes as presented in Table |5| 



For BB5e4 and BBle5 cases in Figure^which are not listed in 
Tables and |^ the plotted points for tev are the actual values 
simulated at fully-refined grid resolution 512 X 1024, adjusted 
slightly to correspond to the converged limit of our results for the 
highest grid resolution of 1024 X 2048, based upon the convergence 
test results in Table [Tl The adjusted values are 0.77 (0.88) times 
the simulation values for resolution 512 X 1024 for cases BB5e4 
(BBle5), respectively. 
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Table 2. Evaporation times and ionizing photon consumption rates per atom for several simulations at highest grid resolution (1024x2048 
fully refined), all for = 9 and rt/(Aa;)njin = 227. 
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Figure 1. Photoevaporation times icv for individual minihaloes 
(top panels) and total ionizing photon consumption per halo atom 
^ (bottom panels) vs. halo mass M (left), redshift z at which I- 
front first encounters the minihalo (middle) and dimensionless 
ionizing flux _Fo (right), for the three different types of source 
spectra, BB5e4 (circles), QSO (triangles) and BBle5 (squares), 
labelled by (M,Fq^z) to indicate the parameters which we kept 
fixed in each case as we varied the quantity on the x-axis. The 
red (green) lines on the two left panels correspond to redshift 
1 + 2 = 15 (1 + 2 = 20). Also, note the different vertical scale of 
the lower right panel only. 



4 RESULTS 

4.1 Evaporation time 

Our simulation results for the evaporation time tov are plot- 
ted for all cases in Figure^ (upper panels). The evaporation 
time increases significantly with increasing halo mass. For 
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Figure 2. Fractional errors of the fitting formulae in equations 
and jSJ compared with the simulation results for each of the cases 
plotted in Figure for the evaporation times (top panels) and 
total ionizing photon consumption rates § (bottom panels), re- 
spectively, vs. halo mass M (left), redshift z at which I-front en- 
counters the minihalo (middle) and dimensionless ionizing flux Fq 
(right), for BB5e4 (blue circles), QSO (red triangles) and BBle5 
(green squares) spectra. 



the smallest minihaloes t^v ~ 10 Myr, much shorter than 
the Hubble time at that epoch, while for the larger mini- 
haloes tav ~ 100 Myr, closer to, but still shorter than the 
current Hubble time. The redshift dependence of the evapo- 
ration time is even stronger than the mass dependence, with 
tcv increasing with redshift. However, due to the relatively 
small range of the relevant redshifts, the overall variation of 
tov with redshift is only by factor of ~ 3 for fixed halo mass 
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Figure 3. Comparison of simulation results for the evaporation 
times and photon consumption rates with predictions of the OTS 
approximation. Ratios of tcv/tsc (top panels) and ^/^OTS (bot- 
tom panels) vs. halo mass M (left), redshift z at which 1-front 
encounters the minihalo (middle) and dimensionless ionizing flux 
Fq (right), for each of the cases plotted in Figure [Tl for BB5e4 
(blue circles), QSO (red triangles) and BBleS (green squares) 
spectra. 



and external flux. Since the Hubble time tn x [1 + z)^^'^, 
the evaporation times of higher-redshift haloes are a much 
larger fraction of the Hubble time, although still generally 
smaller than tn- Finally, the evaporation times depend in- 
versely on the level of the external ionizing flux, becoming 
extremely long {tev ~ few hundreds of Myr) for low flux lev- 
els (Fo < 0.1), and quite short (tov < 30 - 40 Myr) for high 
flux levels {Fo > 10^). 

The logarithmic slope of the dependence of tov on mass 
and the linear slope of its dependence on the initial redshift 
show no signiflcant dependences on the source spectrum for 
the range of masses, fluxes and redshifts we have studied, al- 
though the overall normalization does vary. The dependence 
of tcv on the flux is more complicated, with a slope which 
varies with Fo- 

A simple analytical flt to these results is given by 



tev = AM-^Fg 



C+D login -Fo 



E + F 



Myr 



(5) 



where A = (150,97,128), B = (0.434,0.437,0.465), C = 
(-0.35,-0.357,-0.358), D = (0.05,0.01,0.056), E = 
(0.1,0.3,0.24), and F = (0.9,0.7,0.76) for cases (BB5e4, 
QSO, BBle5), respectively. We have flt the functional de- 
pendences based on all the points plotted in Figure The 
relative errors of these fitting formulae are plotted in Fig- 



ure |5| (upper panels) for all the cases in Figure Table |5| 
shows the comparison between the simulation values and the 
fitting formulae for several of our highest-resolution simula- 
tions. In most cases the errors are <10%, and in all but two 
cases the errors are <20%. 

We have compared our simulation results for 
tev with the OTS e stimat e of this quantity by 
iHaiman. Abel fc Madaul (1200 ll) in Figure ^ The OTS 
approximation, recall, assumes that the sound-crossing time 
of the halo once it is photoheated by ionizing radiation, 
tsc = 2rt/cs, where Cs is the speed of sound at 10'' K, is a 
good approximation to the evaporation time. According to 
equation this means tov.OTS = tsc oc rt oc M^^^ {l + z)~^ , 
so the OTS evaporation time for a fixed halo mass decreases 
with increasing redshift in inverse proportion to the red- 
shift, rather than increasing linearly with redshift, as in our 
numerical results, according to the scaling in equation 
while the increase of tev with halo mass in the OTS 
prediction is somewhat less steep than our result. Finally, 
the OTS approximation completely ignores the significant 
dependence of tev on the magnitude and spectrum of the 
ionizing flux, making tsc a signiflcant underestimate of tev 
(by up to an order of magnitude) for low fluxes and a 
similarly significant overestimate for high fiuxes. 



4.2 Ionizing photon consumption 

The simulation results for the number of ionizing photons 
per minihalo atom required to evaporate a minihalo, 5, are 
plotted in Figure (lower set of panels). The photon con- 
sumption per atom increases steeply with increasing halo 
mass. For the smallest minihaloes, ^ ^ 2, regardless of the 
spectrum of the ionizing source, i.e. each atom on average 
recombines just once during the evaporation of these haloes, 
or, equivalently, is ionized twice by the time it is expelled 
from the halo in the evaporative wind. On the other hand, 
for the larger minihaloes, the photon consumption rate is 
significantly larger, ^ ~ 5 — 8, so each atom in these haloes 
recombines multiple times on average during the evaporation 
of the halo. A BBle5 source evaporates a large minihalo a 
factor of ~ 2 more efficiently in terms of photon consump- 
tion than a QSO or BB5e4 source, although the differences 
in efficiency between different spectra almost disappear for 
the smaller minihaloes. The ionizing photon consumption 
rate grows approximately linearly with redshift, and hence 
is noticeably higher at higher redshifts. Finally, opposite to 
the trend for the evaporation time, the photon consumption 
rate grows with increasing flux, with only 2-3 photons per 
atom needed in the case of low ionizing flux (Fo < 0.1), but 
up to 7-16 needed in the case of high flux levels (Fo > lO'^). 
Therefore, and somewhat counter-intuitively, the low flux 
levels are much more efficient than the high fiux levels, al- 
though the evaporation process takes much longer in the 
former than in the latter cases. 

A good fit to these results is given by 



^ = i + AM^^^^°^^'>'''^F(f'^'^^°'^^'>^° 



F + G 



'1 + 
. 10 . 



,(6) 
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where A = (4.4,4.0,2.4), B = (0.334,0.364,0.338), C = 
(0.023,0.033,0.032), 

D = (0.199,0.24,0.219), E = (-0.042,-0.021,-0.036), 
F = (0,0,0.1), G = (1,1,0.9) for the BB5e4, QSO, and 
BBle5 spectra, respectively. The relative errors of these scal- 
ing laws are plotted in Figure |21 lower panels. In most cases 
the errors are <10%, and in all cases but one the errors are 
below 20%. 

We have compared our simulation results for ionizing 
photon consumption with the predictions of the OTS ap- 
proximation in Figure 01 The OTS prediction for ^ IS given 
by equation (|1J with / = 1. The number of recombina- 
tions per atom, ^ots — 1, scales with mass and redshift as 
M^/^(l -I- zf. While the logarithmic slope of this OTS de- 
pendence on the halo mass is roughly similar to that of our 
simulations according to equation © , the increase of ^ with 
redshift of the OTS prediction is much stronger than that 
of the simulation results. In addition, the dependence of ^ 
on the ionizing flux and spectrum is quite strong in the nu- 
merical results, while it is completely ignored in the OTS 
approximation. Unlike the evaporation time, t^v, however, 
the overall normalization of ^ in the OTS approximation is 
overestimated by at least an order of magnitude in all cases, 
as compared with the simulation results, and by more than 
two orders of magnitude in some CEises, particularly at high 
redshift and/or low ionizing flux levels. 

4.3 Evolution of the neutral mass fraction 

As described in Paper I, the first phase of the encounter be- 
tween an intergalactic I-front and a minihalo is the weak, 
R-type phase, during which the I-front enters the halo and 
photoionizes an outer layer of halo gas, predominantly on 
the side facing the ionizing source, before the gas has time 
to respond hydrodynamically and move. This phase ends 
when the I-front slows to the R-critical speed and makes a 
transition to D-type, thereby trapping the I-front until it 
can burn through the remaining neutral, shielded minihalo 
gas and evaporate it. We identify the fraction of the origi- 
nal minihalo gas mass which was photoionized in this initial 
R-type phase when the front became R-critical as Ma- We 
have plotted in Figure|l]the values of Ma for our simulation 
cases in Figure □ In Paper I, we derived an analytical ap- 
proximation which can be used to estimate Mcr, called the 
Inverse Stromgren Layer (ISL) approximation, where ISL is 
the ionized boundary layer of the minihalo defined by the 
"outside-in" Stromgren length for each impact parameter 
(see Paper I for details). In the static, ionization-equilibrium 
approximation, the mass contained in this boundary layer, 
MisL, which is readily obtained for our TIS model by numer- 
ical integration as shown in Paper I, can be used to estimate 
Mcr as the complement of Misl, Mct = A/tot — AfisL- 

In Paper I, we compared our simulation results for the 
ionized mass fraction at the transition from R-type to D- 
type for a lO^M© halo exposed at = 9 to a source of flux 
_Fo = 1 and found an excellent agreement with Mish- Here 
we extend the comparison to determine how well that agree- 
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Figure 4. Neutral mass fraction of minihalo at the moment of 
I-front trapping, Mcr. Points show results from simulations with 
BB5e4 (circles), QSO (triangles) and BBleS (squares) spectra, re- 
spectively. For comparison, we plot the analytical prediction for 
this quantity derived from the Inverse Stromgren Sphere (ISL) 
approximation in Paper I (with argument on x-axis shifted as de- 
scribed in text), in which Mcr^analyt 

= 1-AfisL/A^tot (solid line), 
and a convenient fit to the ISL approximation and simulation re- 
sults, according to equation (dashed line). 

ment holds up when we vary the halo parameters, source flux 
level, and redshift. We find that the dependence of the simu- 
lation results on these parameters follows that predicted by 
the ISL approximation fairly well, where the latter predicts 
that AfisL should depend upon (Af, Fo, z) only in the com- 
bination FoM~'^''-^{l + z)~"\ We find, however, that the fuU 
range of simulation results for Mcr is better fit by this ISL 
scaling law if we multiply the quantity FqM'^^'^ z)^"' in 
the original ISL approximation by 10. The excellent agree- 
ment between this modified ISL scaling law for Mcr and the 
simulation results is shown in Figure 01 We find that the fol- 
lowing analytical expression provides a reasonable fit both 
to the ISL curve (with shifted argument) and the simulation 
results, 

Afcr.approx = 0.9954 exp(-0.0013a:'' ) , (7) 

where a: = 5 + Xog^oWoM-^''' {1 + z)-^]. 

In Fig. we show the evolution of the neutral mass 
fraction, M„/Mcr vs. t/tcv as the minihalo evaporates, for 
all of the cases plotted in Figure^ This time- variation of the 
neutral mass fraction of the evaporating minihaloes is ap- 
parently a nearly universal function, when the neutral mass 
M„ it) is expressed in units of Afcr and t is expressed in units 
of tcv for each case, independent of the halo parameters and 
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Figure 5. Evolution of the neutral mass fraction Mn/Mcr from 
the end of the initial R-type phase of the I-front and its conver- 
sion from R-type to D-type to the final evaporation of the halo, 
for simulations with BB5e4 (top), QSO (middle) and BBle5 (bot- 
tom) spectra (solid lines). We also show the best fits in each case, 
as described in the text (long-dashed red lines). 

external flux, but allows for some dependence on the spec- 
trum of the external source. The scatter in all cases is largely 
due to the decrease of the flux level in time due to cosmo- 
logical expansion, for a source located at flxed comoving 
distance from the target minihalo. Accordingly, the relation 
is tightest in the QSO case, since in these runs the evapo- 
ration times are shorter. The average shape of the neutral 
mass time-dependence function is well-fit by 

M„{t) = M,, (l - , (8) 

where A = (1.07,1.03,1.03) and B = (2.5,2,2) for the 
BB5e4, QSO and BBleS spectra, respectively (red dashed 
lines in Fig. O. 



5 IMPLICATIONS FOR REIONIZATION AND 
CONCLUSIONS 

The widespread presence of minihaloes filled with neutral 
gas during the epoch of reionization made encounters be- 
tween the weak, R-type intergalactic I-fronts which ionized 
the IGM and these minihaloes a common occurrence. This 
interaction affected both the minihaloes and the I-fronts pro- 
foundly. The trapping of the I-fronts by minihaloes which 
cover the sky as seen by the more massive haloes which 
are generally believed to have been the primary sources of 



cosmic reionization slowed the advance of those I-fronts in 
the IGM, at large, and wasted some of the ionizing radi- 
ation which would otherwise have been available to ionize 
that IGM. The minihaloes, in turn, were transformed by the 
photoevaporation which followed their trapping of the inter- 
galactic I-fronts, into barren, dark-matter haloes devoid of 
baryons. In order to explore this process and its implications 
for cosmic reionization, we performed a large set of high res- 
olution gas dynamical simulations which include radiative 
transfer, the first in their kind. In Paper I, we discussed our 
results in detail for the illustrative case of a 1O''M0 halo 
which is overrun by an intergalactic I-front at Zi =9, from 
a source of flux _Fo = 1, for different source spectra. Here we 
have summarized a much larger set of simulations designed 
to quantify the dependence of minihalo photoevaporation on 
the halo mass, source flux level and spectrum, and redshift. 
In Paper I, we discussed the comparison between our simu- 
lation results for IO'^Mq haloes and previous related work. 
Here, we extend this comparison to include the broader set 
of cases presented in this paper. 

We noted in Paper I that the approximate analysis of 
iBarkana fc Loeb (^1999) which was used to argue that pho- 
toevaporation of minihaloes affected a signiflcant part of 
the baryon fraction collapsed into haloes before the end of 
reionization was inconsistent with the resu lts of our detailed 
numerical simulations for IO^A^q haloes. iBarkana fc Loebl 
|1999) had modelled this process by a static approximation 
like the ISL approximation in Paper I, to determine what 
portion of the halo's mass is shielded from external ionizing 
photons and what portion is exposed to photoheating. They 
then assumed that only that portion which was instanta- 
neously heated in this way, enough to unbind it from the halo 
gravitational potential well, would evaporate from the halo. 
From this static approximation, they concluded that only 
minihaloes smaller than ~ few x 10^ — 10® Mq (depending 
on the redshift and the source spectrum) were evaporated 
completely during reionization, while the larger minihaloes 
would have retained a significant fraction of their gas (up to 
~ 50% of the gas for haloes of mass ~fewx lO'^Af©). In Paper 
I we demonstrated that this model does not account prop- 
erly for the dynamics of photoevaporation, and, as a result, 
it fails to anticipate the fact that even minihaloes as large in 
mass as M = 10'^A/q must ultimately evaporate completely. 
The larger set of simulations we have presented here show 
further that this neglect of dynamics by Bark ana fc Loebl 
( 1999) in estimating the dependence of photoevaporation on 
halo mass is not correct; all minihaloes exposed to ionizing 
radiation during reionization evaporate completely. All of 
the gas initially in a minihalo would eventually have become 
unbound, outfiowing with speeds of Uwind = 20 — 40kms~^, 
leaving behind just a dark halo. The static approximation 
fails because it ignores the dynamical nature of photoevap- 
oration. As evaporation proceeds, layers of gas are continu- 
ously stripped away, exposing the gas layers within to the 
ionizing radiation. Thus, although the inner halo region self- 
shields initially, all the gas is eventually "unshielded" and 
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photoheated to T > 10* K, well above the halo virial tem- 
perature, and, hence, boils out of the halo. 

In Paper I we pointed out that simulations of the pho- 
toevaporation of 10^ haloes by an intergalactic 1- front 
during reionization did not support the suggestion by ICenI 
(|2001j) that external ionizing radiation can cause the implo- 
sion of the minihalo gas, leading to globular cluster forma- 
tion. We did not observe such an effect in our simulations 
reported there. The broad range of cases considered here 
allow us to extend this comparison to the full range of pa- 
rameter space expected for minihaloes exposed to the effects 
of externally-driven I-fronts during reionization. In no case 
ha ve we yet fo und evidence of the implosion effect predicted 
bvlCenI poo ill . 

Our simulations show that the number of ioniz- 
ing photons per minihalo atom needed to photoevapo- 
rate the minihalo is typically in the range ^ ~ 2 — 
10, significantly larger than the value of ^ ~ 1 pre- 
viously estimated_t^^ IGM 
jGnedinI l2000al: iMiralda-Escude. Haehnelt fc Reed l200Cl) . 
On the other hand, wc have also shown that recent claims 
iHai man. Abel & Madau 2001) that the photoevaporation 
of minihaloes can require up to hundreds of ionizing pho- 
tons per atom are very significantly overestimating ^. 

We have quantified here the dependence of evapora- 
tion times and photon consumption rates on the halo mass, 
source flux level and spectrum, and redshift of encounter 
between I-front and minihalo. We have provided convenient 
fitting formulae for our simulation results as functions of 
the parameters (M, Fo,z), which should be useful in further 
analyses of the role of this process in cosmic reionization. 
To go beyond this step, it is necessary to determine the dis- 
tribution in time and space of these parameters. The effect 
on reionization will depend upon the luminosity function of 
source haloes and the mass function of minihaloes as they 
evolve in time and fluctuate in space, including the clus- 
tering of both the sources and minihaloes. These may de- 
pend, in turn, upon the feedback of reionization on galaxy 
formation, so it is likely to be necessary to solve the prob- 
lem of reionization in some detail before the true impact of 
minihaloes as screens and photon sinks is known. That is 
well beyond the scope of this paper. Here, instead, we shall 
provide a simple first estimate of the average effect of mini- 
halo photon consumption on the photon budget required for 
reionization, based on the mean Press-Schechter (PS) distri- 
bution of minihaloes, as follows. 

Let the ionizing photon consumption rate per minihalo 
atom be given for a single minihalo by ^(M, z, Fq). Then its 
average over all minihaloes at redshift z is given by 



Ptot/coll,MH 



" i[M,z,F,)M'-!^dMM 



where dn{M,z)/dM is the well-known PS mass function of 
haloes. 
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Figure 6. Collapsed fraction in minihaloes /coll,MH (top panel) 
for the ACDM universe, mean photon consumption per minihalo 
atom, ^, due to photoevaporation of minihaloes (middle panel), 
and mean excess photon consumption per total atom (i.e. all 
atoms, including both minihaloes and IGM) compared to the 
nominal requirement of one per atom (bottom panel). Results 
are for ionizing fluxes _Fo = 0.1 (long-dashed), 1 (solid), 10 (short- 
dashed) and 100 (dotted), for case BB5e4. 



is the collapsed mass fraction in minihaloes, ptot is the 
mean matter density at the corresponding redshift, Mmin 
is the Jeans mass in the uncoUapsed IGM, Mj, and Mmax ~ 
M{W^K), is the halo mass corresponding to virial tempera- 
ture lO** K. Assuming that ^(Af, z, Fo) is given by the fitting 
formula in equation @, we can calculate ^{z,Fo) for any 
given redshift z and ionizing flux Fo- Results for Fo = 0.1, 
1, 10 and 100 and 7 ^ 1 + z ^ 21 are shown in Figure]^ 
for case BB5e4. We see that ^ depends strongly on the ion- 



izing flux (^ oc Fq 



0.199-0.042 logio Fo 



in equation 0), but is 



(10) 



approximately independent of the redshift. Thus, the con- 
tribution of minihalo photoevaporation to the mean global 
photon consumption rate is largely dictated by the ionizing 
flux levels and the current collapsed fraction in minihaloes 
and can require up to ~ 1 additional ionizing photon per 
atom in the universe (both minihalo and IGM atoms) to 
finish reionization, as compared with the requirement when 
minihaloes are neglected, for the flux levels we considered 
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here. Hence, minihalo photoevaporation can potentially dou- 
ble the number of ionizing photons required to reionize the 
universe. 

We defer the complete treatment of minihalo and source 
bicis which would account for the spatially and temporally- 
varying fluxes during reionization to a future paper (Iliev, 
Scannapieco & Shapiro, in preparation). A brief sum- 
m ary of our first results along these lines can be found 
in llliev. Scannapieco. Shapiro fc R.agal l)2004h . where it is 
shown that the statistical bias which causes the minihalo 
mass function to be enhanced around source haloes relative 
to the mean PS mass function boosts the relative contribu- 
tion of minihaloes as photon sinks and compensates for the 
declining mean collapsed fraction at high redshift. 

This shows that minihalo photoevaporation will be an 
important feature of reionization, causing a modest, but sig- 
nificant slow down of the global I-fronts and delaying their 
overlap. Only a more detailed, self-consistent treatment of 
the global process for a given reionization scenario will deter- 
mine if reionization was photoevaporation- ciommaied, since 
the influence of minihaloes is strongly dependent on both the 
flux and, through the evolving collapsed fraction, the red- 
shift at which the I-fronts encountered the minihaloes. Ad- 
ditionally, minihalo formation in some places and at some 
epochs might be partially suppressed by e.g. an early X- 
ray background, double reionization, or reduced small scale 
power. Such different scenarios might be distinguished obser- 
vationally by detecting the fluctuations of the re dshifted 21- 
cm emission or ab sorption lines from minihaloes (llliev et alJ 
I2OOI l2003l : ICenI 12003') or the kind of absorption line sig- 
natures at UV and optical wavelengths discussed in Paper 
I. 
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